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We have studied the measurement accuracy of the forward-backward asymmetry of the 
l^^ ' e'^e" — >■ ti process in the 6-jet mode at -^5 = 500 GeV with the ILD detector. In the 

^SJ , analysis the vertex charges of 6-jets were used to distinguish t from i in each event. 

The distribution of the cosine of the reconstructed polar angle of so identified t or t 
showed a clear forward-backward asymmetry. After the correction for charge misidenti- 
fication the forward-backward asymmetry was determined to be A%g — 0.334 ± 0.0079 



X 



O 



for 500 fb with the beam polarization combination of P(e , e ) = (-1-30%, —80% 



Mh' demonstrating a very good statistical accuracy (~ 2%) even in the 6-jet mode. 



1 Introduction 



►v^ , The forward-backward asymmetry of the e^e^ — > ti process is sensitive to the tiZ and tij 

l/-\ ' couplings and serves as a probe for new physics, which may appear as anomaly in these 

ff^ , couplings. 

^^ ■ Needless to say we have to distinguish t from i in each event in order to measure the 

\l^ , forward-backward asymmetry, A^g . In the lepton-f4-jet mode, it is straightforward because 

lO ' the lepton charge tells the charge of the leptonically decayed W, and hence identifies its 

^^ , parent to be either t or i. In the 6-jet mode, however, we need some other way to separate 

t and i, which is non-trivial. Nevertheless it is worth investigating the feasibility of A^g 

measurement in the 6-jet mode, since it has a major branching fraction of 46%. In addition, 

the kinematical fit works better in the 6-jet mode than in the lepton-|-4-jet mode, where 

►v> , a large energy is taken away by the neutrino from the leptonically decayed W . The 6-jet 

^ • mode might, hence, be advantageous when the influence of beamstrahlung is significant. 

5^ , In this paper, we use the vertex charges of b jets to distinguish t from i in each event. 

The measurement of the vertex charges requires a high performance detector system as well 

as a sophisticated reconstruction algorithm. The forward-backward asymmetry in the 6-jet 

mode is therefore a very good measure for the overall performance of a detector system, 

hence being included as one of the ILC Lol benchmarks [1] . The benchmark conditions are 

the center of mass energy of -y/s = 500 GeV, an integrated luminosity of 500 fb~^, and a 

beam polarization combination of P(e+, e^) = (-1-30%, —80%). 

The ILD [2] is equipped with an unprecedentedly excellent vertex detector, which allows 
efficient 6-jet charge identification with the vertex charge. The vertex charges of b jets were 
reconstructed by the LCFIVertex algorithm [3]. 

This paper is organized as follows. After we briefly describe our analysis framework in 
section [21 we move on to the vertex charge reconstruction and show how well we can identify 
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each jet as 6 or 6 in section[31 We then apply this to t/t identification for the determination of 
the production angle distribution {dN /d cos 9t) in sectional After the correction for charge 
misidentification we derive the Ay,^ and discuss the result in comparison with the Monte 
Carlo truth in section [S] Section |6] summarizes our analysis and concludes this paper. 



2 Analysis framework 

In general, Monte Carlo (MC) simulation consists of the following steps: event generation, 
detector simulation, event reconstruction, and data analysis. 

All of the MC event samples, both the signal and the backgrounds, used in this study 
were produced in the StdHep [S] format by a SLAC team jjj as common inputs to Lol 
studies, using WHIZARD 1.4 [B] for generating parton 4- momenta and PYTHIA 6.2 [7] 
for parton-showering and hadronization. The beam energy spread and the beamstrahlung 
were properly taken into account using the spectrum generated with Guinea-Pig [5] for the 
default ILC design parameters at ^/s = 500 GeV. 

The final-state particles output in the StdHep format from the event generation step 
were passed to a Gcant4-based detector simulator called Mokka [9] and swum through the 
ILD detector to create exact hits in trackers and calorimeters. 

These exact hits were smeared or digitized, if necessary, depending on the detector 
elements in the first part of MarlinReco [10] [11]. The pattern recognition was done for 
the smeared tracker hits, separately in the TPC and the silicon detectors and so found 
track segments were then linked together and fed into a Kalman-filter-based track fitter in 
the second part of MarlinReco. From the fitted tracks and the calorimeter hits, individual 
particles were reconstructed as particle flow objects (PFO) with a sophisticated particle flow 
algorithm called PandoraPFA [T2] in the third part of MarlinReco. 

These PFOs were forced to cluster into 6 jets for the signal and all background events 
with Durham jet clustering algorithm [T3] in the fourth part of MarlinReco. 

The next step is heavy flavour tagging with LCFIVertex [3]. The LCFIVertex package 
consists of two parts. The flrst part is to search for secondary and tertiary vertices inside 
each jet and locate them, thereby determining the decay length, the transverse mass, and the 
momentum at each of these vertices. In the second part these quantities are translated into 
the impact parameter joint probability and the highest two impact parameter significances, 
which are used as inputs into a neural net (NN) trained with jet samples having and 1 or 
more secondary vertices. Each reconstructed jet is then assigned with the three NN outputs, 
corresponding to &-, c-, and 6c-tags. 

Once a bottom-fiavoured jet is identified we can determine whether it is 6 or 6 by mea- 
suring the vertex charge. We will discuss this in detail in section [3] 

The reconstruction of e^e~ — > tt events in 6-jet final states has been studied extensively 
in the context of the top quark mass measurement. The top quark mass in the 6-jet mode 
is one of the benchmark observables for ILC Lois and the ILD study has been reported in 
reference [2]. Similar to the top quark mass measurement, the measurement of the forward- 
backward asymmetry requires a correct jet-parton association. In this paper we employed 
the same reconstruction algorithm as those used for the top quark mass measurement. For 
the reconstruction, therefore, we refers the readers to the above reference. 
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3 Vertex charge reconstruction and its performance for a single 
6-jet 

In principle, we can tell which three jets are from top and the rest from anti-top in the 6-jet 
final state by identifying either a 6/6-jet or a c/c-jct from W~^ /W~ decay as shown in Fig.[TJ 
In this analysis, however, we used only 6/6-jcts for the top/anti-top separation. 




Figure 1: Schematic diagram showing fully hadronic decay channel of e"'"e — > tt process 

A bottom quark hadronizes into a B-hadron, which flies over a finite distance thanks 
to its large ct, making a secondary vertex significantly away from the primary vertex and 
hence identifiable by the vertex detector. We define the vertex charge as the sum of the 
charges of the charged tracks associated with the secondary vertex. If the charged tracks 
are reconstructed and associated perfectly, the vertex charge is equal to the charge of the 
primary i?-hadron, from which the charge sign of the b/b is uniquely determined. 

In practice, however, there is no perfect vcrtexing, and the resultant distribution of the 
vertex charges of charged B-hadrons will have a finite width and hence their charges might 
sometimes be mis-identified. By the same token, the charges of neutral B-hadrons might be 
mis-identified, causing confusion in bottom charge sign identification. 

Figure [2] shows the distribution of the MC level vertex charge, which is defined by the 
charge of a i?-meson {B^ or i?°) involved in the jet tagged as a 6-jet. 

Inspection of the figures tells us that we can separate B+s from B~s by selecting jets 
with negative or positive vertex charges with some contamination from B'^s. 

Figure [3] shows the distribution of the reconstructed vertex charges of the jets which are 
^-tagged. In this distribution the 6-quark (6-quark) contribution is shown by hatched blue 
(red). 

The charge sign identification efficiency for a single 6-jet using the vertex charge is 28% 
with a purity of 75%. Notice that only 40% of the b-jets hadronize into charged B-hadrons 
and hence maximum efficiency one can hope for is 40% in this method. 
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Figure 2: Reconstructed vertex charge distributions for (a) charged B-mesons and (b) 
neutral i?- mesons as separated using MC truth information (PDG particle ID code). 
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Figure 3: Reconstructed vertex charge distributions for 6-quarks (hatched red) and 6-quarks 
(hatched blue), as identified using MC truth information, shown together with that of all 
the reconstructed 6-jets (solid gray). 
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4 t and t identification and determination of the production angle 
distribution {dN / d cos 6 1) 

Let's call two reconstructed top systems ti and i2 and 6-tagged jets associated to them 61 
and 62 , respectively, The identification of t and t is performed by using the vertex charges 
of &i and 62 as follows. We define Ci {i = \ and 2) as the vertex charge of hi and C = ci — C2 
as the event charge. If C is 0, the event is thrown away as we cannot tell t from t. If C is 
positive, ii is t, while if C is negative, ii is t. The typical distribution of the event charge 
C is shown in Fig. U) 
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Figure 4: Distribution of event charge: C = ci — C2. 

Figures [S]-(a) through -(i) show the distributions of the production angle of ti for 9 
possible combinations of the signs of the vertex charges: ci and C2, for the samples where 
both W bosons decayed into light quarks. 

We can see a clear forward-backward asymmetry in the case of C 7^ (Figs. [S]-(a) to 
-(f)). Notice also that in the first row (-(a) to -(c)) ti is t while in the second row (-(d) to 
-(f)) ti is t, and hence showing opposite asymmetries. On the other hand, in the case of 
C = 0, for which the charge signs of both 6-jets were undetermined or they were inconsistent 
(both 61 and &2 had the same sign), there is no forward-backward asymmetry visible. In 
each figure the contribution from the events with wrong charge sign is hatched red, showing 
an opposite forward-backward asymmetry. Hatched blue is that from the combinatorial 
background in which the reconstructed 6-jet candidates were not actually 6-jets, showing no 
forward-backward asymmetry. 

The combinatorial background depends on the flavour into which W bosons decay. Fig- 
ures [6]-(a) through -(i) are the same figures as Figs. [5]-(a) through -(i), but plotted for the 
samples where one of the two W bosons decayed into a c/c-quark. Figures [T]- (a) through 
-(i) are similar figures for the samples where both W bosons decayed into c/c-quarks. We 
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can see clearly that the combinatorial background grows with the number of c-jets in the 
final states, since the probability of mis- identifying charm as bottom increases. 

If C is positive, ti is t, while if C is negative, ti is t. If we can assume that <i and ^2 are 
back-to-back, the production angle of t is obtained from the angle of ti as 

cos^f = (Tc ■ I cos 6*4 J, (4.1) 

where ac is the sign of C . In order to test this assumption we compared the distributions 
for ii with those of ^2 after reflection (hatched green) in Figs. [S] |6l and[7]-(d) to -(f). This 
comparison confirmed the assumptiorO, allowing us to use Eq. ()4.ip to combine all the figures 
with C 7^ 0. The selection efficiency of this cut is 69%. 

The differential cross section for e^e~ — ^ tt and consequently its forward-backward 
asymmetry [A^pg) depend on the center-of-mass energy of the ti system, which, in turn, 
depends on the amount of initial state radiation and beamstrahlung. In order to make A^pg 
well-defined, we hence rejected events with vs significantly less than 500 GeV, by requiring 
the acoUinearity between ti and <2 to be less than 8°. This final cut discarded 24% of the 
events so far survived. The overall selection efficiency is 20% for fully hadronic tt events. 

The resultant production angle distribution is shown in Fig. [H Of the final sample 83% 
have correctly identified signs of top quark charge. 

5 Determination of A'^p^ 

The production angle distribution Fig. [8] is distorted by the charge mis-identification. The 
distortion can be corrected by using the following formulae: 

dNohsiO) = p{0) ■ r]{9) ■ C ■ da{9) + P{tt - 6) ■ r]{7r - 9) ■ C ■ da{-K - 0) + dB{e) 
iNohsiTT -9) = p(7r - 9) ■ i]{tt - 9) ■ C ■ daiir -9)+ p{e) ■ rj{e) ■ C ■ da{0) + dB{iT ~ 9) 

(5-2) 
where rj{9) is the acceptance, p{9) {p{9)) is the probability of correctly (wrongly) assigning 
the charge sign at production angle 9, and C is the integrated luminosity. 
p{9) and p{9) are plotted in Fig. [HI as the black and red lines, respectively. 
The figure shows no ^-dependency, allowing us to set p{9) = p and p{9) = p. Solving for 
da{9), we thus obtain the following formula for the differential cross section: 

p ■ jdNMO) - dB{9)) - p ■ jdNMn -9)^ dBjir - 9)) 

'^"^'^ = {p^-p-)m-c ■ ^'-'^ 

Figure [TO] shows the production angle distribution after the correction for the charge 
mis-identification. 

The acceptance function rj{9) is in general 0-dependent. In our case, however, it turned 
out at that the 6'-dependence was negligiblq^. The acceptance hence cancels out in the 
calculation of the forward-backward asymmetry, resulting in 

^^ k<e<.n''^i(^)^L,2<e<.dN{9y ^'^ ^ 



* Strictly speaking, Eq. 14.11 1 does not hold on an event-by-event basis because of initial state radiation 
and beamstrahlung. On average, however, ii and t2 are back-to-back, allowing us to merge the six cases 
with C 7^ 0. 

' Notice that for the reconstructed t or t, being a 3-jet system, there is no acceptance hole even at 6 = 

or TT. 
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Figure 5: Distributions of the reconstructed polar angle of the identified top quark in fuUy- 
hadronic tt events, which are tagged as (Qh^Qha) == i^) (+'^); (b) (+,0), (c) (0,—), (d) 
(-,+), (e) (-,0), (f) (0,+), (g) (0,0), (h) (+,+), and (i) (-,-) for the sample in which 
both W bosons decayed into light-quark pairs {bbuddu sample). 
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Figure 6: Distributions of the reconstructed polar angle of the identified top quark in fuUy- 
hadronic tt events, which are tagged as (Qh^Qha) ~ i^) (+'^); (b) (+iO), (c) (0,—), (d) 
(-,+), (e) (-,0), (f) (0,+), (g) (0,0), (h) (+,+), and (i) (-,-) for the samples in which 
one of the two W bosons decayed into a c/c-quark {bbudsc/ bbcsdu samples). 
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Figure 7: Distributions of the reconstructed polar angle of the identified top quark in fuUy- 
hadronic tt events, which are tagged as (Qh^Qha) == i^) (+'^); (b) (+,0), (c) (0,—), (d) 
(-,+), (e) (-,0), (f) (0,+), (g) (0,0), (h) (+,+), and (i) (-,-) for the sample in which 
both W bosons decayed into c/c-quarks {bbcssc sample). 
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Figure 8: Distribution of the reconstructed polar angle of the identified top quark in fuUy- 
hadronic ti events. The contributions from events with the wrong charge (red) and the case 
where the 6-quark is mis-identified are also shown (blue). 
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Figure 9: cos 9ti dependence of the charge-sign identification (black) and mis-identification 
(red) probabilities. 
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Figure 10: cos^ti distribution after the correction for the charge mis-identification. 



where dN{9) = r]-C-da{e). 

From Fig. [10] we finally obtain 



A^PB = 0-334 ± 0.0079 (stat.), 



(5.5) 



where the center of mass energy of ^/s 



and a beam polarization combination of P(e+, e 



500 GcV, an integrated luminosity of 500 fb 

(+30%, -80%). 



6 Summary 

We have studied the measurement accuracy of the forward-backward asymmetry with the 
ILD detector for the e+e^ — > ii process in the 6-jet mode at -y/s = 500 GcV. In the analysis 
the vertex charges of 6-jets were used to identify t and i. The efficiency to tag the vertex 
charge of a 6-jet was 28% with the purity of 75%. Having two 6-jets in each event, the 
probability to identify t and t in the event was 71% with the probability of wrong charge 
assignment of 12%. The measured angular distribution was corrected for wrong t/i charge 
assignments. From the number of events in forward and backward hemispheres after the 
correction, we obtained A^pg = 0.334 ± 0.0079, where quoted error is statistical only. 
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